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Abstract 



Oh ' Motivated by the simulation of stable random fields, we consider the issue of 

^ ■ discrete approximations of independently scattered stable noise. Two approaches are 

', proposed: grid approximations available when the underlying space is and shot 

noise approximations available on more general spaces. Limit theorems stating the 
convergence of discrete random noises to stable white noise are proved. These results 
are then applied to study moving average spatial random fields with heavy-tailed 
I innovations and related limit theorems. A second application deals with discrete 

^ ■ approximation for Brownian Levy motion on the sphere or on the euclidean space. 

i/^ . Key words: stable white noise, stable fractional noise, linear filtering. 

o : 1 Motivations 

a^ 

^ ' Stable integration is a basic tool in the theory of stable random process. Indeed, if W 

denotes a stable random noise on the measurable space {E,S) with control measure 
^ ■ m and {ft)t&T a kernel such that ft G L°'{E,£,m) for any t ^ T, then the random 

process defined by 

Xt= [ ft{x)W{dx), t E T 

is a stable process, and the path properties of the process can be deduced from the 
properties of the kernel (see for instance [IB] chapter 10). For example, the first ex- 
amples of self-similar stationary increments (SSSI) stable process were constructed in 
such a way: using a moving average kernel, Taqqu & Wolpert [19] and Maejima p3] 
constructed the linear fractional stable motion; the harmonisable fractional stable 
motion was proposed by Cambanis & Maejima [2] as the stable integral of the har- 
monisable kernel. More recently, Cohen and Samorodnitsky [4] propose a new class 
of SSSI stable process defined as the stable integral of a random kernel, the random 
kernel being the local time of a fractional Brownian motion. Generalizing the notion 
of selfsimilarity for random field in higher dimension, operator scaling stable random 
fields were defined and constructed by Bierme k, Scheffler [1] as stable integrals of 
high-dimensional kernel with suitable scaling properties. 
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For the purpose of simulation of these processes, we need a general theory of dis- 
crete approximation of stable random measures. An approach based on Lepage's 
series also called shot noise series was developped: generalized shot noise series were 
introduced for simulation in [15], further developments were done in [16] and p2] and 
a general framework was developed in [3j. In this paper, we propose a new approach 
for the discrete approximation of stable random noise. We propose and developp two 
alternatives: in the case when is a stable white noise on W^, a grid approximation 
can be used (see Theorem 13.1(1 : for more general spaces, an approximation of the 
stable white noise by a Poisson random measure is discussed (see Theorem 13. 2p . 
We apply our results on grid approximation of stable random noise on M!^ to the study 
of stable noise obtained by linear filtering from i.i.d. random fields in the domain of 
attraction of stable distributions. This extends the results of Kokoszka & Taqqu 
[U [HI Uni Hi] for one dimensional sequences such as ARMA or FARIMA with stable 
inovations to spatial random fields. We prove a limit theorem for the suitably rescaled 
random noise (see Theorem 14. ip . Two cases occur: if the coefficients of the random 
filter decrease quickly, the limit random field is stable white noise and the dependence 
vanishes in the limit; if the coefficients decrease slowly like a power function, the limit 
random field is fractional stable random noise with long range dependence. 

Another application of the Poisson approximation for stable white noise is given 
in the framework of stable Levy motion on the sphere S'^ or on the euclidean space 
W^. These processes mimick the simple covariance structure of standard Brownian 
motion on R to more general metric spaces (see tl2j)- As a direct application of our 
results, we give a Donsker's type Theorem for stable Levy motion on the sphere or 
on the euclidean space. 

The paper is organised as follows. In section 2, we remind the reader of general 
results on random noise, stable white noise, Poisson random measures and their con- 
vergence. In section 3, the convergence of discrete random measure to stable white 
noise is proven for two different models: grid approximations and Poisson approxi- 
mation. In section 4, two applications are exposed: linear filtering of i.i.d. random 
field and Donsker's type theorem for Levy stable motion on the sphere or on the 
euclidean space. Section 5 is devoted to the proof of our results. Technical results on 
convergence of deterministic functions are gathered in an appendix. 

2 Stable random noises and Poisson random mea- 
sures 

A random noise on {E,£) is a generalized random field {W[f])f^jr indexed by a 
linear subspace T of the space of real- valued measurable functions on E, verifying 
the linearity property: for all ai, a2 G K and /i, /2 G ^- 

W[aih + a2f2] = aiW[h] + a2W[f2\ a.s. (1) 
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Using the Cramer- Wold's device, this imphes that the finite-dimensional distributions 
of the process (W^[/])/e:r are determined by its one-dimensional distributions. The 
set of functions T is thought as the set of integrands. To emphasize the analogy with 
usual integration, we use the notation 

W\f\ = ! f{x)W{dx). 
Je 

\i 1a ^ ^ for some measurable set ^ C the value W^(^) = VF[1a] is thought as 
the measure of A. However, each realisation of W does not define an usual signed 
measure because it needs not to have finite variations. That's why some author 
(including myself) prefer the term random noise to random measure, this last term 
being reserved for measure-valued random variables. 

We recall in this section the definition and properties of stable white noises (in- 
cluding Gaussian white noises) and Poisson random measures. 



2.1 Independently scattered a-stable random noises 

Let {E,£,m) be a measured space with m a ci-finite measure. For a > 0, we de- 
note by = L"{E,£,m) the space of measurable functions /:£'—> M such that 

|/|°£im < oo. For / G L°, let = (/^ |/|°dm)^/" . If a < 1, this defines a 

norm and is a Banach case. This is no longer the case if < a < 1. 

The stable distribution of index a £ (0, 2] and parameters a > (scale), G [—1)1] 
(skewness) is denoted by 5c(c7, z^). For the sake of simplicity, we will always assume 
that = if a = 1. The distribution Sa{cr, v) is defined by its Fourier transform 

A(^) = exp(-cj°|0|" (l-izye(0)tan^)) , 0gM (2) 

where e(a) = if a > 0, e{a) = — 1 if a < and e(0) = 0. 

The a-stable random noise Wa on {E,£) with control measure m and skewness 
function v : {E,£) [—1, 1] is the stable random process defined on some probability 
space and indexed by L", verifying the linearity property [T] and with one 

dimensional marginal distributions given by 

where 



cjj = (^J |/|"dm^ ^ and afuj = J 



Mf)\f\"dm. (3) 



Evaluating the process Wa at functions {fi)i<i<k with pairwise disjoint supports 
yields independent random variables {Wa[fi])i<i<k' we say that the random noise 
Ma is independently scattered or white to qualify this property. 

In the Gaussian case a = 2, the parameter z/ is irrelevant and we retrieve the usual 
Wiener integral, which is an isometry from Lp'{E) onto some Gaussian subspace of 
L^(r2,.F, P). In the case 1 < a < 2, the stable integral induces an isometry from 
L°^{E) onto some subspace of stable random variables on {Q,T,¥) equipped with the 
covariation norm (see 11811. 
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2.2 Poisson random measures 

Poisson random measures are measure-valued random variables and they induce a 
random noise on their sets of integrands. In the following presentation, we focus on 
the point of view of random noises, a general reference is |7j. 

Let {E,S,n) be a measured space with n a cr-finite measure. We denote by Sq 
the class of measurable sets with finite n-measure. A random measure N on (E, £) 
is called a Poisson random measure with intensity n if and only if for any finite 
collection (^i)i<i<n of pairwise disjoint sets in the random variables (-/V(^i))i<i<n 
are independent with N[Ai) following a Poisson distribution with mean n{Ai). The {a 
priori random) space L^{E,£, N) of integrable functions with respect to A'' is almost 
surely equal to 

rN = {f ■E^R \ \f\MeL\E,n)] 

and N induces a random noise on !Fn. We note N[f] = fdN for / S Tn- This 
random noise is characterized by the linearity property ([T|) and the one dimensional 
marginals given by their characteristic functions 

E [exp(i6'Af[/])] = exp (e^^-^ - l)dnj , f eJ^N, OeR. 

2.3 Convergence of random measures 

In this paper, we are mainly concerned with convergence of random noises {fJ-h)h>o to 
stable random noise W, either independently scattered (Theorem I3.ip or fractional 
(Theorem 14. We shall now give a precise definition of convergence of random noises. 

For /i > 0, let fih be a random noise with integrand space .F^^ and let W be some 
random noise with integrand space J-w- Let furthermore ^ be a subspace included 
in all the J^^^^s, h> and in J^w- 

Definition 2.1 We say that fih converges weakly to W on as h ^ and write 

Hh =^ W as /i ^ 0, 

if the finite dimensional marginal distributions of {fih[f])fej^ weakly converge to those 
ofmf])fe^- 

Using the linearity property of the random noises and the Cramer- Wold's device, 
we easily see that convergence of finite dimensional distributions is equivalent to con- 
vergence of one dimensional distributions, which is in turn equivalent to convergence 
of characteristic functions. Although straightforward, this will be of constant use and 
is stated for future reference in the following proposition: 

Proposition 2.1 The following statements are equivalent, where convergence is meant 
as h ^ 0.- 

1. converges weakly to W on T , 

2. for any f ^ T , /U/i[/] converges weakly to W[h], 
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3. for any f eJ^ and 6" G R, E [exp(z6'/x,,[/])] converge to E [exp(z6'l^[/])]. 



To illustrate this definition, we give a first example of weak convergence of random 
noises. 

Example : Let {E, d) be some metric space and £ the borehan (T-algebra. Let 
m and (mh)h>o be finite measures on E and v and be measurable functions E 
[—1,1]. Define W (resp. Wh) as the a-stable random measure with control measure 
m (resp. rrih) and skewness function v (resp. Vh)- Let Ci,{E) denotes the space of 
bounded continuous functions on E. Note that from the assumption that m is finite, 
Cb{E) CLTw = L"{E, £, m), and similarly Cb{E) C = L"{E, £, nih) for all /i > 0. 

Then the following two statements are equivalent, where convergences are meant 
as ^ — > 0: 

1. the random noises Wh converge weakly to W on Cb{E), 

2. the control measures nih converge weakly to m and the skewness measures 
Ufidrnfi converge weakly to udm. 

The proof is direct once we recall that the a-stable distributions Sa{cyh^^h) converge 
to Sa{(y-, I') (with a > 0)if and only if {ah, I'h) {o; 

3 Discrete approximations of white stable noises 

We introduce two methods for approximating an independently scattered stable ran- 
dom noise on {E,£). The first one, available when the underlying space is = R'^ 
and the control measure is the Lebesgue measure, relies on a grid approximation 
hZ" C W when the span h of the grid goes to zero. The second one, avalaible on a 
general space E, relies on the approximation of the stable random noise by suitable 
Poisson random measures or shot noise, when the intensity A goes to infinity. 

3.1 Grid approximation of white stable noise on 

In this section, E = and Wa denotes the independently scattered a-stable random 
noise on M.'^ with Lebesgue intensity and constant skewness v. We propose a discrete 
approximation of Wa based on the grid approximation /iZ*^ C M*^ with span h > 0. 

The construction uses a family ^ = {^fe, k G Z*^} of real random variables satisfying 
the following assumption: 

Ha the ^fe's are i.i.d. random variables in the normal domain of attraction of the 

stable distribution »Sq((T, v). 

We suppose furthermore = if a = 1. Recall that the normal domain of attraction 
of the stable distribution (Sq,(o", z/, 0) consists in the random variables Y such that 

n 

n~^^" ^ 5a (cr, z/) as n — oo 

i=l 
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where the Yi are i.i.d. random variables distributed as Y. If a = 2, H2 holds if and 
only if the ^fc's are i.i.d. random variables with ^[Ck] = and IE[^^] = 2(7^. The stable 
distribution 52(0", 0) is then the normal distribution with mean and variance 2(7^ 
(the skewness parameter v is irrelevant in this case). In the case < a < 2, from the 
known characterization of the domain of attraction of stable distributions (see [6]), 
the ^fc's must satisfy the following tail assumptions as x — > +00: 

P(efc >x)^ px~" and P(Cfc < -x) ~ qx-", (4) 

where the constants a, i^,p,q satisfy 

POO 

{p + q) / t""sintdt, 

/>oo 

ip-q) / t""(l - cost)dt. 
JO 

Note that equation ([4]) is equivalent to the fact that ^ — r belongs to the normal 
domain of attraction of the a-stable distribution Sa{cy,v) for some r G M, but we 
assume here furthermore that r = 0. 



a 



(j"utaxi{-Ka/2) 



We propose the discrete approximation of Wa obtained by replacing Wa{dx) by 
lhik'^h{k+i'i){x)dx on each cell h{k + I*^), k G . Here I*^ = [0,1)"^ and 7/1 = 

(T^^/i'- Q More formally, let = ^^h{0 be the random measure on R"^ absolutely 
continuous with respect to Lebesgue measure with random density 



dx 



{x) =lh^ ^kih{k+i'i)ix). (5) 



The random signed measure fih defined by ^ induces a random noise on the 
integrands set 



•^Mh ^ {-^ ^ ^loci^'^) ; j^^ fdfJ'h converges a.s. | . 



(6) 



Note that in this definition, only the semi-convergence of the integral J^d fd^h or 
equivalently of the series 

Ih'y^ik I f{x)dx 
Jhik+It) 



is required almost surely. 

The following scaling relation is worth noting: 

/i,[/(c.)]=c-"Vh[/(.)], c>0 (7) 

whenever these quantities are well-defined. 

Our first result precise the integrands sets We need to introduce the space 
P"(M'^) of locally integrable function / G Ll^iM'^) such that /(x) = o(|x|-'') for 
some r] > a~^d. 
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Proposition 3.1 Let < a < 2 and suppose that ^ satisfies assufnption Hq,. 
For h > 0, consider the random measure f^h defined by equation Then, 

1. //!<«< 2, L-(M'^) C^^,. 

2. 7/0 < a < 1, P"(M^) c T^^. 

We then consider weak convergence of the random measures {fJ^h)h>o- To unify the 
notation, we introduce J^" = L°(R'^) if 1 < a < 2, and J^" = P"(M°') if < a < 1. 

Theorem 3.1 Let < a < 2 and suppose that ^ satisfies assumption Hq. 
Then as h —> 0, converge weakly to on . 

Remark: In the case 1 < a < 2, we obtain convergence on L"(]R'^) the full natural 
set of integrands for Wa- In the case < a < 1, we have to restrict on the smaller 
subspace P"(M'^) because L° is not included in (see the remark after Lemma [52] 
in the Appendix.) 

Remark: We prove in fact a slightly stronger result 

^ as /i ^ if A ^ / in To.. 

See Proposition 15. II for a precise statement including the definition of convergence in 
P"(M'^) if < a < 1. This diagonal convergence can be used to prove the following 
interesting corollary. 

Let C(M'') denote the space of continuous functions on M*^. 

Corollary 3.1 Let < a < 2 and suppose that ^ satisfies assumption Hq,. 
Then the random signed measure fl^ on defined by 

converges weakly to Wa on T"' n C(]R'^) as h —>■ 0. 

These results can be used to propose simulation of a-stable processes represented 
as stable integrals. Another approach developed in [3] is to use the representation 
of stable integral in Lepage's series. For the purpose of simulation, it is important 
to provide error bounds or speed of convergence for the method. We now consider 
the speed of the convergence of the random noise fih to the stable integrals Wa 
in Theorem 13.11 Better rates of convergence are expected if the random variables 
^ satisfy Hq with distribution Sa{cr,i^) rather than in the domain of attraction of 
the stable distribution. In this case, the discrete approximation fj,h[f] has the same 
distribution as 

fih[f] = I iM{x)Wa{dx), 

where iphf is the discrete approximation of / obtained by replacing / by its mean 
value on each cell h{k+ [0, k (see equation [TQ] below). In this representation 
corresponds to ah~'^/°'Wa[h{k+ [0, !['')] ~ 5a((T, u). As a consequence, the random 
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variable can be represented on the same space and we can focus on LP convergence 
and give error bounds. Furthermore, during the simulation, the kernel / must 
sometimes be replaced by some approximated kernel /m- For example, when / has 
infinite support, /m might be the restriction of / on a bounded domain. Or /m 
might be a piecewise constant approximation of / so that the mean value of /m on 
h{k + [0, reduces to fMihk). 



Proposition 3.2 Let < a < 2 and /m € be an approximation of the kernel 
/ G JT" and let 



and Ca,p is the absolute p-th moment of the distribution u). 
Ln the case a = 2, the result holds for all p > 

Remark: Note that iphfM fM in as ^ — ^ and explicit error bounds can be 
given under suitable regularity condition on /m- 

3.2 Shot noise approximation of white stable noise on a 
general space 

In this section, (E,S,m) is a general measured space with m a o"-finite measure and 
Wa denotes the independently scattered a-stable random noise on E with control 
measure m and constant skewness v. We propose a discrete approximation of Wa 
based on shot noises, i.e. on Poisson random measures. 

Let G be some distribution on R belonging to the normal domain of attraction 
of the stable distribution Sa{a,iy). For A > 0, let Nx{de,d^) be the Poisson random 
measure on x M with intensity measure nx{de,d£,) = Xm{de)G{d^). It induces a 
random noise on E by the heuristic formula 




Then, for every < p < a 




where 7^ = cr 
that 



a . More formally, let J^^^ be the set of functions f : E ^M. such 




(8) 



We consider the random noise fix on J^^^ defined by 




(9) 



Our result is the following: 
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Theorem 3.2 Suppose G belongs to the normal domain of attraction of the stable 
distribution Saia,i^). Then, for any A > 0, L"'{E,£,m) C J^^^ and as X ^ oo, nx 
converge weakly to Wa on L'^{E,£,m). 

Remark: In the case when m is finite and normahzed to m{E) = 1, we can 
replace the Poisson random measure A'^;^ by the binomial random measure Nn = 
a"~^n~a ^"^^ (5(^. where {Ci,ei)i>i is a sequence of i.i.d. random variables with 
distribution G{d^)m{de). We then consider the random noise 

n 

f^nif] = o-"^n"^ ^iifiei). 

1=1 

Similar results as those of Theorem 13.21 hold: for any n > 1, L°'{E,£,m) C J^^^ and 
jln weakly converge to Wa on L°'{E,£,m) as n — > oo. The proof is very similar to 
that of Theorem 13.21 and will be omitted. 



4 Applications 

We propose two applications of the above results about convergence to stable random 
noise. The first application based on Theorem 13.11 is an analysis of random noises 
arrising from linear filtering of i.i.d. random fields. Interessant fields appear in the 
limit that we call fractional stable random noise. The second application based on 
Theorem 13.21 proposes a Donsker's type theorem for stable Levy field on the sphere 
or on the euclidean space. 

4.1 Spatial linear filtering 

In this section, we consider a random measure fi^ of the form 

=^hYl ikih{k+i'i)i^), (10) 

where % > is a normalisation constant and the field ^ = (Cfc)^^^^ is not any more 
i.i.d. but stationary with spatial dependence (possibly long range dependence). More 
precisely, we consider the case when the field ^ is obtained from a random field ^ 
satisfying by linear filtering. The linear filter is given by c = (cA;)^g^d G r and 
the random field ^ by 

L = {c*c)k = ^ ck-iCi- 

The corresponding assumption is 

Ho ^ = ^ * c with ^ satisfying Hq and c G 

Here = ^"(Z*^) is the space of sequence u = Wk)^^'£d such that X^fcgz'' l^^l" < 
and * denotes the convolution product of sequences. We will see that if ^ satisfies Hq, 
and c G then * c is defined almost surely. 
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The integrands set associated to the random measure jlh is 



J^jx^ = |/ G LlM'^) ; j^^ fdP'h converges a.s. | . (11) 

We consider the convergence of the random measure jlh to some a-stable noise. 
According to the behavior of c at infinity, the hmit field can be either a white stable 
noise or a fractional stable noise (to be defined below). In the former case, spatial 
dependence disappear in the limit and weak dependence holds, whereas in the latter 
case, spatial dependence persists in the limit and strong dependence holds. For the 
sake of simplicity, we focus on the case 1 < a < 2. As in section 13.11 Wq, denotes a 
stable white noise on M"^ with Lebesgue control measure and constant skewness v. 

Theorem 4.1 Let 1 < a < 2 and suppose that ^ satisfies assumption Hq,. 
For h > 0, consider jlh the random measure given by equation / flOj) . 

1. Suppose c £ i-^ and let % = a~^hS'^ ^-i)d (j _ (^^^^^ c^). 

Then for any /i > 0, C !F^^ and as h ^ 0, fih converge weakly to CWa on 

2. Suppose that c has regular variations at infinity of order —(3 G {—d, — a^^d) in 
the sense that 



lim sup 



t^cyt^]-p{x) =0, (12) 



where p : {0} ^ M is some locally integrable homogeneous function of order 
-(3 (hyp ??). Let% = a-^h'^~^^-P. 

Then for any h > 0, H CI JT^^ and as h ^ 0, ^h converge weakly to Wa,p 
on n L^, where Wa,p is the fractional stable random noise defined by 

Wa,p[f]=Wa[f*p] , /GL^nL" 

We now give some interessant properties of the random noise Wa,p- 
Proposition 4.1 Let 1 < a < 2 and p be as in Theorem \^.l\ 

1. The random noise Wa,p on n L^ is a-stable, stationary and a~^d — (5 + d 
self-similar. 

2. In the case a = 2, W2-p is the centered Gaussian random noise with covariance 
function 

E {W^,p[fWaM) = L ™, fi^)Kix, y)giy)dxdy, f,geL'nL^ 

where K : the kernel given by 

K{x,y)= p{x - z)p{y - z)dz. 
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Remark: The notions of stationarity and selfsimilarity are invariance properties 
of the random noise under the transformation group of translations and dilatations 
respectively. The formal definition is given in equation ?? and ??. 

Remark Self-similarity can be stated using the notion of renormalisation groups. 
Let he a linear space of functions on M*^ closed under dilatations, i.e. for any h > 
and / € jr, f{h.) £ JF. Let u G M. For /i > let T^"'' be the transformation acting 
on a random noise on by 

{Tj:'^^W)[fi.)] = h-W[f{h.)] , fGT. 

The group relations T^^^ = Id and T^"^ o T^"^ = T^^, holds and the group of trans- 
formation (r^"'*)/i>o is called the normalisation group of index u. Our results have 
a nice interpretation in terms of normalisation group. Consider first Theorem 13. 1[ 
The scaling relation jTl implies that the random measure fih is equal to Tjf '^^ Hi 

and the Theorem states that the random noise T^" '^^ f^i converges to Wa as /i ^ 0. 
Furthermore, Wq, is a fixed point of the renormalisation group of index a~^d, i.e. 
jX" f^)^^ _ meaning that is (a~^d)-selfsimilar. In the same way, con- 
sider Theorem 14.11 part 2. Here jlh = Tj:° ^'^'^^ jli and the Theorem states that 
^'^'^^ fii converges to Wa,p as /i — > 0. Furthermore, Wa,p is a fixed point of the 

renormalisation group of index a~^d — f3 + d, i.e. T^" ^~^'^^Wa,p = Wa,p, meaning 
that Wa,p is Q!~^(i — /3 + d-selfsimilar. 

4.2 Convergence to Levy Brownian motion on the sphere 
or on the eucUdean space 

The Brownian motion parametrized by a general metric space was introduced by 
P.Levy. It is defined as follows: let O be a fixed point in a metric space (M, d), a 
Brownian motion parametrized with the metric space (M, d) and with origin O is a 
centered Gaussian process {B{m))m£M such that : 

- B{0) = almost surely, 

- B{m) — B{m') has variance d{m,m'). 

There does not always exist a Brownian motion on (M, d) for an arbitrary metric 
space (M, d) but constructions of such Levy Brownian motion have been proposed 
in the case when M is a sphere, a euclidean space or an hyperbolic space (see [12J). 
These constructions are based on Gaussian white noise and have been extended to 
stable white noise in order to get stable processes with interesting properties. 

First consider the case when M = S"? is the unit sphere in R'^^^ and d is the 
geodesic distance. Let O be a fixed point on the sphere 5"^, for example the northern 
pole and let W2 be a Gaussian white noise on the sphere with control measure ds the 
normalized uniform measure on the sphere 5"^. For m £ S"^, let Hm be the hemisphere 
centered at m defined by 

= {m' E 5" I d{m,m') < ^}. 
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Then the process {B{m))m&si defined by 

B{m) = (HoAH^) , m £ 

is a Brownian motion on {S'^,d). Here A denotes the symmetric difference. More 
generally, the Levy stable motion on the sphere S"^ is defined by the same formula 
with the Gaussian white noise W2 replaced by a symmetric cc-stable random measure 
Wa with control measure ds. 

When the metric space {M,d) is the euclidean space R'', a Brownian motion is 
constructed in the following way (Levy-Chenstov construction). Let W2 be a white 
Gaussian noise on 5^ x with control measure ds x dr. A pair (s,r) represents 
the hyperplane in with equation < x, s >= r and -S^ x is thought as the set 
of hyperplanes that do not contain the origin. For m G define Vm as the set of 
hyperplanes that separate O and m, i.e. 

Vm = G S« X ; < r «s,m>}. 

Then the process {B{m))^^i^<i defined by 

B{m) = W2 {Vm) , meW 

is a Brownian motion on Replacing the Gaussian white noise W2 by the 

symmetric a-stable random measure Wa with control measure dsdr, we obtain the 
symmetric a-stable Levy-Chenstov random field, which is a a~^-self-similar stationary 
increments process. 

We apply our results on convergence of Poisson random measures to stable noises 
to obtain a Donsker's theorem for stable Levy motion on the sphere S'^ or on W. 

Theorem 4.2 Suppose G is in the normal domain of attraction of the distribution 
(Sa((T, 0) for some a > 0. The following convergence hold in the sense of finite dimen- 
sional distributions as A 00. 

1. Let Nx{ds, d^) be a Poisson random measure on S'^xM. with intensity XdsG{d^). 
The random process B\ on S'^ defined by 

BaM = 0Fa-iA-V« /■ ^lHoAHMNxids,d^) , meS'i 

weakly converges to the symmetric a-stable Levy motion on the sphere 5'*. 

2. Let N\{ds,dr,d^) be a Poisson random measure on 5"^ x x R with intensity 
XdsdrG{d£,). The random process Bx on defined by 

Bx{m) = a-^-V" f ^lv^{s,r)Nx{ds,dr,dO, m G R« 

JSixR^xR 

weakly converges to the symmetric a-stable Levy motion on the euclidean space 
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Remark : The random fields Bx can be easily simulated, which is not clear for the 
limit field B. For example in the case of Levy motion on the sphere S"^, simulating the 
Poisson integral can be made as follows. Draw T according to a Poisson distribution 
with mean A, conditionally to T = t, draw (sj, '^j)i<i<t identically distributed with 
distribution dsG{d(,) and let 

t 

i=l 

Note that the quantity x S Hm is easily calculated since x S Hm if and only if 
< x,m >> 0. In the case of Levy motion on the sphere M'^, the method is just the 
same once we observe that simulating B\{m),m G in a bounded domain D requires 
the knowledge of the Poisson random measure on a bounded domain of S'^ x and 
hence only of a finite number of random points. 



5 Proofs 

5.1 Preliminaries on stable distributions 

We recall some known facts about normal domains of attraction of stable distribution 
(see [6]). Let ^ belong to the normal domain of attraction of the stable distribution 
Sa{(y^v). Then, the following estimate holds for its characteristic function as — > 

\{e) = E = \{e) + od^l"). (13) 

where A is given by Furthermore, in the case < a < 2, the tail estimate ([4]) 
implies that there exists C > such that for any s > 

Var[a{|^|<4] < C^s^-" and E[|e|l{|5|<4] < Cs^-". (14) 

5.2 Proof of Proposition 13.11 

The random measure defined by ^ is linked with the discretization of the space M'^ 
by /iZ*^, /i > 0. Introduce the operator 



Define the random signed measure u on l!^ by 



It defines a random noise on the integrands set J-'i, defined by 

= |/ : Z'^ ^ M; sumj^^j^dfkik converges a.s. | 
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where once again convergence means semi-convergence of the integral. We have the 
formal relation 

^^h[f] = ihhMM] 

from which we deduce that / G .F^^ if and only if iphf G J^^- But we will see below 
that = (see Lemma O) and that 'iph{L°') C ^ if 1 < a < 2 and C 
if < a < 1 (see Lemma [5?2l in the Appendix). This proves Proposition 13.11 □ 

Lemma 5.1 Suppose a G (0,2] and ^ satisfies Hq. Then, Fy = 

Proof of Lemma 15.11 

This is a direct application of Kolmogorov's three series Theorem (see [6j). The 
random series Ylike.Z'^ fkCk with independent summands converges almost surely if 
and only if for any s > 0, the following three numeric series converge: 

{\fA\<s}] < 



and 

J2 E [fk^kl{\M,\<s}] converges. 

Note that the convergence of the first series for small s implies that — > as 
/c ^ oo, so we do suppose fk^Oask^oo. In the case a = 2, this implies 

Var [fkCkl{\M„\<s}] ~fc^oo 

and hence the second series converges if and only if / G Suppose now / G 
Cebycev's inequality implies 

p[iMfei>5]</|E[ei]s-2 

so that the first series converges. Applying Borel-Cantelli's Lemma, this in turn entails 
that the event limsup{|/fc^A;| > s} has probability 0. Hence, the set {k; \fk£,k\ > s} is 
almost surely finite, so that in the third series almost all terms vanish and the third 
series converges. 

In the case < a < 2, we deduce from equation ([4]) that 

>s]~(p+'7)iAr5-" (16) 

and hence the first serie converges if and only if f G Suppose now that f £ i". 
Equation p4l) implies the second and third series converge. □ 



5.3 Proof of Theorem 13.11 and of Proposition 13.21 

Proof of Theorem [37T] 

Using Proposition 12.11 it is enough to prove the convergence of one dimensional dis- 
tribution: for / G .F", 

^^h[fHWa[f] as h^O. 
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We will prove in fact a stronger result that will be useful in the sequel and we consider 
diagonal convergence. Recall that T'^ = L'^ if 1 < a < 2 and J-°^ = I?" if < a < 1. 
If a G [1,2], is a Banach space when endowed with the norm ||.||l" and the notion 
of convergence //^ — > / in is clear. In the case < a < 1, we need the following 
definition of convergence in P". 

Definition 5.1 We say that {fh)h>o converge to f in "D" if the following two condi- 
tions hold: 

• for any compact K C , fn^K converges to fix in , 

• there is some rj > a^^d such that fh{x) = o(|x|~'') and f{x) = o(|x|~^) as 
X ^ oo uniformly in h > 0. 

Theorem 13.11 is then a direct consequence of the following Proposition. 

Proposition 5.1 Suppose ^ satisfies Hq and let {fh)h>o converge to f in . Then, 
the following diagonal weak convergence holds: 

l^hifhHWM] as h^O. 

We prove convergence of the characteristic functions. Let 6 £ M.. From the 
definition of stable white noise, 

E [exp{ieWa[f])] = exp {-af\er{l - iu}e{e) tan(7rQ/2))) (17) 
with CTj and Vf given by ([3|). On the other hand, we can rewrite 

lihlfh] = a-^h'^''" iMh){hk), (18) 



where V"/! : L\^^ — > L\^^ is the linear functional defined by 



{'il^hf){x) = h-'' / f{u)du, (19) 

and [x\h = h[h~^x\. Note that 'ij^hf is the approximation of / when the space M"^ is 
discretized by hTJ^ and the function / is replaced by its mean value on each cell of the 
form /i(A: + /"'), k E Z'^. With these notations, the characteristic function of ^Jih[fh\ is 
given by 

E [exp(i0/x^[A])] = n ^ {cT-^h''~'\^hfh){hk)e) . (20) 
Here we have used equation (fTSll . the independence of the ^^'s and the almost sure 



convergence of the random serie (|18l) implying the convergence of the above infinite 
product. 

First step: we begin to show that as ^ 0, 

E [eMiOl^hUh])] = n ^ (a-^h''''\Mh){hk)e) + oil). (21) 
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To see this, we estimate the difference and use the following inequality : let {zi)i^i 
and {z[)i^i two families of complex numbers in D such that the products Y\i£i 
Yliel convergent, then 



Using equation (f2Q]) . this yields 



feez 



< Yl I AK(A:))-AK(A:)) 



(22) 



with 

Equation ([T3|) implies that the function g defined by g{0) = and 

giv) = \v\-'' \\iv) - ~X{v)\ , v^O, 
is continuous and bounded so that for any k G Z'^, 

|A {uh{k)) - A {uh{k))\ = g{uh{k))\uh{kT. 
In order to obtain an uniform estimation, define the function g : [0, +oo) 

by 

g{u) = sup \g{v)\. 

\v\<u 



[0,+oo) 



Note that g is continuous, bounded and vanishes at 0, and that for any A; G Z*^ such 
that \uh{k)\ < e, 

I A K(A;)) - A K(A;)) | < g{e)\u„{k)r. (23) 
Let e > 0. Equations ([20]), ([22]) and ^ together yield 

E [exp{i9iih[fh])] - n ^ [a-'h''-'''ii,hfh)ihk)9 

< 5(e) Yl M\uh{k)\<e} + 2 ^{\uUk)\>e}- 

Here we have used equation ([23]) to bound \X{uh{k)) — X{uh{k))\ from above whenever 
|^i/i(^)| ^ e, and the trivial upper bound 2 otherwise. Now we remark that 

E Wh{krh\u,ik)\<e} < E Kikr = ^''^mi^hfkn^ 
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is bounded since from Corollary 15.11 \tphfh\°' |/|° in whenever //i — > / in 
By the continuity of g in 0, g{e) is small for e small enough. Furthermore, 



^ '^{\uH{k)\>e} - ^ j^^^{\i,hfh{x)\>ea\e\-^h-'--'d}'^^ 



and this quantity goes to as /i — > because Corollary 15.11 implies that \tphfh\'^ — > |/|" 
in wheneyer fh^f'm. !F°' and hence the familiy /i > is uniformly inte- 

grable. These estimates imply equation (pTI) . 

Second step: we proye that as ^ — > 0, 

J] \(a-^h''~"\^hfh){hk)d) =exp(-^7^|0r(l-iz.;£(0)tan(7ra/2))) +o(l). 

(24) 

To see this, we take here advantage of the exponential form of A and write the l.h.s. 
of (El as 



<Q> 



= exp (^-l^r \{Mh){x)\^dx + ii/tan(7ra/2)r-> y |(^;,A)(x)|<">dx ) (25) 

where u^"^^ = e(n)|n|" denotes the signed power function. 
Corollary 15.11 implies that 



\{^hh){xTdx= / |/(x)rdx + o(l) 



and 



l(V^hA)(x)|<">dx= / |/(x)|<">dx + o(l). 
This together with equation ([25]l implies equation (f24l) . 



Conclusion: The convergence of the characteristic functions of iJih\fh\ given by ([20]) 
to the characteristic function of W[f] given by (flTl) is a direct consequence of equa- 
tions ([2T]) and (f24l) . This proves the weak convergence Hh{fh\ =^ ^[/]- Proposition 
15.11 and Theorem 13.11 are proved. □ 

Proof of Corollary [3^1 

Note that formally jlh[f] = where fh is the function defined by fh{x) = /([x]/i). 

The result follows from Proposition 15.11 and from the fact that if / G H C(M'^), 
then h^fin T'' . □ 
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Proof of Proposition 13.21 

The proof is rather straightforward once we observe that p'h[fM]—Wa[f] = Wa[(phfM— 
f] has the same distribution as \\(phfM — fWia^a where Xq, has distribution 5q,(1, i^). 
Hence 

E[|Ah[/A/] - = ||^,/m - f\\lE[\X^n 

and = E[|X„|P]. □ 

5.4 Proof of Theorem [X2 

We first prove that for any A > 0, the integrand set contains L°'{E,£,m). Let 
/ : £^ — > M be such that J^\f\"dm < oo. We have to prove that equation ^ is 
satisfied and compute 

/ \^f{e)\Alnx{dC,de) = A / (|/(e)|E[|e|l|/(e)5|<i] + P[|/(e)C| > 1]) m(de) 

< A / C|/(e)|"m((ie) < oo 



JE 

where the last inequahty is a consequence of equations (fTBl) and (fT4l) . This proves the 
required inclusion. 

We now prove the convergence of the random measures to the stable random 
noise. Following Proposition 12.11 we compute the characteristic function of Nx[f] for 
f eL^{E,£,my. 

kJexr / 

where ^(u) = e*" — 1. Performing integration with respect to G, the integral in the 
right hand side rewrites 



E 



{ejxCf{e)))Xm{de)G{dO, 
with "^du) = J^"^{u^)G{d^). We now prove the following asymptotic result: 
hm / ^G{Olxf{e))\m{de) = - S |0/(e)r (1 - zve(0/(e)) tan(7ra/2)). 

A— >oo 



E JE 

This is a consequence of Lebesgue's convergence Theorem. Indeed 7a — > as A ^ oo 
and the asymptotic behaviour of at given by equation (fT3]) yields 

*G(^7A/(e)) = -|^/(e)r (1 - MOf{e)) tan(7ra/2)) + o(l). 

Furthermore, there is some C > such that < C|0|" for some C > and this 

implies the domination condition 

where the r.h.s. is integrable with respect to m{de) and does not depend on A. Hence 
we have proved that 



lim E [eicp{ieNx[f])] = exp - / |e/(e)|"(l - ve{e f{e))taii{T:a/2))m{de 
and this is precisely the characteristic function of Wa [/] • □ 
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5.5 Proof of Theorem 14.11 

In view of Lemma ISTTl the conditions c G i'^ and ^ satisfies Hq ensure that = is 
well defined. This convolution relation implies that the random measures fj,h and flh 
defined by equations ^ and (fTOjl respectively are linked by a simple relation. This 
is the object of the following Lemma. In order to unify the two cases considered in 
Theorem 14.11 we denote the integrands set by J^a with J^a = in the first case and 
= L° n in the second case. We define c G by Ck = C-k, A; G Z'^. 



Lemme Under the assumptions of Theorem 14. II (both cases), J^" C J^^^ and the 
relation flh[f] = tJ'h['>Phf] holds for any / G T'^ with 

i^hf = -M{i^hf)*c), (26) 
7/1 



where 



and 



(l>h: * , . . (2^ 



Proof: 

We have formally, 



k&Z'^ 

= 7/1 XI {'^hf)kCk-lil 

{k,i)eZ'^xZ'^ 

= —iJ'h[(t>h{{'4^hf)*c)\. 
Ih 



Furthermore, this formal computations are valid because the assumptions of Theorem 
14.11 entail (V'/t/) * c G and hence ip'f^f G C JF^. Recall indeed that Holder 
inequality implies that u * f G ^" as soon as u G £^ and u G Here in the first case 
(resp. in the second case) c G and V^/i/ £ (resp. c G and il)hf G l^). □ 
Proof of Theorem I4.lt 

Theorem 14.11 is a direct consequence of the above Lemma, of Theorem 13.11 (or rather 
of Proposition ISTTl about diagonal convergence) and of the convergence of in 
as /i — > 0. Lemma [531 in the appendix states indeed that il^f^f Cf converges to Cf 
(resp. to / *p) in the first case (resp. second case). □ 
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Appendix: some results on space functions 

The following Lemma gathers some useful properties of the linear operators tph,''Ph 
and (ph defined by equations (fT9]) . ([271) and ([28]) respectively. It is standard material 
for a > 1 but some extra care is needed when < a < 1. 

Lemma 5.2 1. The relations iph = 4'h ° ''Ph, i^h^ <Ph = Id jd and iph ° "iph = i^h 
hold. 

2. (phin C i° and for any u £ UhuU- = h'^/°\\u\\ec . 

3. If a > 1, %ph{L'^) = (-'^ and iphi^'^) C L"^ and iph induces a linear projection on 
L". Furthermore, if fh ^ f in L" as h ^ 0, then iph{fh) f in L"- 

4. If < a < 1, iph{Va) C r, and i'hi'I^a) C (Va). Furthermore, if fh ^ f 
strongly in T>a as h ^ 0, then tphifh) f strongly in V^. 



Remarks: Extra care is needed in the case < a < 1 because iphiLjoc ^ -^") is not 
included in as the following example shows: for 7 > and 5 € M consider the 
function / G Lj^^ defined by f{x) = if x E [k,k + k~'^] for k > 1, and f{x) = 
otherwise. Then it is easily seen that / G if and only if — 7 < — 1, and that 
V'l/ G L°' if and only if a{6 — 7) < —1. Comparing these to conditions, we can even 
find / G n such that Vi/ ^ That's why we need some stronger condition 
on / in the case < a < 1 and we introduce the space of rapidly decaying functions 



Proof of Lemma 15. 2t 



1. This is routine verification using the definitions. 

2. Let uei°. Then, 

/ \4'huix)\°'dx = / Imh-iJ^'dx 

JbbsRd JbbsR'' 

and this implies (phU S L"' and = ^'^''"ll^*!^"- 

3. Let a > 1 and f £ L". Using Jensen inequality with the convex function 



y ^ 12/1, 



< 



h{k+I'i) 



f{u)du 
\f{uTdu 
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This implies that V'h/ e and \\iphf\\i'- < /^"''/"II/IIl" • Together with point 1. 
and 2., we easily prove that iph induces a linear projection on L". 
We now prove that V'ft converge weakly as a linear operator to the idendity. First 
consider the case when / is continuous with compact support. Then from the 
continuity of /, for every x G Md {iphf){x) f{x) as /i ^ 0. The assumption 
that / is compactly supported implies that Lebesgue's convergence theorem 
applies and iphf — > / in L". 

The space of continuous compactly supported function is dense in and the iph 
are (uniformly) bounded linear operator. Hence the convergence tphf — > / holds 
for every f E L°'. Diagonal convergence is a consequence of the equicontinuity 
of the iph which all have linear norm less than 1. 

4. Let / G P". Prom the definition of P", there is some integer M > and 
C > and some tj > a~^d such that if |x|oo > M then \f{x)\ < C\x\'i^ . As a 
consequence, if /i < 1 and |a;|oo > M + /i, we have h{[h-^x] + I"^) C B{0,MY 
and 



\i^hf{x)\' 



h-'^ / /(n)dn 
h{[h-^x]+I<i) 



This implies j^d \ilJhf{x)Y^dx < oo since ar] > d and hence iphf G L°'. 
Now let fh ^ f strongly in V^. Prom the above discussion and the definition 
of strong convergence in T>", we see that a similar inequality holds uniformly in 
h < 1: there is some M > 0, C > and rj > a~^d such that if |x|oo > M then 
|/(x)| < C|x|oo' and also for any h < 1 \ ffi{x)\ < C|x|oo''. As a consequence, for 
any e > 0, there is some M > such that for any h < 1 

[ li^hfhixTdx < e (29) 

On the other hand, for fixed M > 0, we prove that \tphfh\°''^B{o,M) uniformly 
integrable. Let K = B{0,M + 1). Then fhlx ^ /Ix in L\ and point 3. 
implies that il^hifh'^K) — > f^K in L^. Since the restrictions of iphifh^K) and 
ijhfh to B{0,M) are equal, this imphes i'4'hfh)'i-B{o,M) converge to flB(OM) 
L^. This in turn implies that |'i/'/i//i|"l_B(o,M) is bounded in for p = > 1, 
and hence uniformly integrable. This together with equation ([29]) implies the 
uniform integrability of the sequence IV'n/nl"- D 

Corollary 5.1 Let a G (0,2]. Suppose {fh)h>o converge to f in T"^ as h ^ 0. Then 
\MfhT ^ l/r inL^ ash^O and also {Mfh))<''> ^ m 

Proof of Corollary I5.lt 

We consider the convergence of {'>phfh{x))^'^^ to /^"^ in L^, the other case being 
treated in the same way. Pirst of all, as above, for any e > there is some M > 
such that for any h < 1 



\x\oa>M 



(VhA(:E))<">-/(x)<"> dx<e. (30) 
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For fixed M > 0, the application L^{B{0,M)) R, f ^ /ij(o,Af) is 
continuous. Hence the convergence {tphfh)'^B{o,M) ~^ f^B{o,M) ™ imphes 



n— ►oo 



lim / (V';.A)(x)<">dx= / (^;,A)(x)<">dx. 
'b{o,m) Jb(o,m) 



Together with equation ([30I1 . this yields the result. □ 



We now consider the convergence of tpf^f defined by (|26|1 . Lemme Under the 
assumptions of Theorem 14.11 

1. In the first case, for any / E L", Tpf^f — > Cf in as /i ^ 0. 

2. In the second case, for any f ^ V- r\ , ipf^f — > / * p in as /i — > 0. 
Proof of Lemma 15. 5t 

1. Consider first the case when c £ We easily see that 

i''hf= E cM{- + hk). 

For each fixed k, the sequence of function i^hfi' + hk) converge to / in (use 
the fact that according to Lemma 15.21 point 3, iphf ^ / in and that the 
translation operator / — > /(. + hk) is continuous and converge to identity as 
h^O). Then for any fixed M, 

X] (^k^hf{- + hk) ^ Ck)f 

\k\<M \k\<M 

in as /i — > 0. Furthermore, since c G the remainder 

II ^M{- + hk)\\L'^ < I Y icfci I ii/iil« 

|A:|>M \|fc|>A/ / 

is small when M is large. Thanks to these estimates, we prove that ipf^f — > Cf 
in L". 



2. Consider the second case when c is such that 



lim sup 



t— >oo 



x\=l 



0, 



and note that ipf^^f can be rewritten as 



rhf{x)= ^hf{[x\h-y)Ph{y)Ay = {{i^hf)*Ph){[x]h) (31) 
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with 



PhW = h —C_[h-ly] 



K 



Using equation [121 we see that for any y / 0, ph{y) Piu) as /i — > 0. 
More precisely, equation [T2l entails that ph^B —>■ Ph^B in and p^ls^ — > Ph^B'^ 
in L", where -B denotes the unit ball -6(0, 1) in M*^ and B'^ its complementary 
set. Using Lemma [5^ and the assumption / G n L", we have also iphf — > / 
in and in L°. The bilinear application x L°' ^ L°', (51,52) i—*^ 51 * 92 is 
continuous and 

*52|Ilc < ||5i|Ili \\92\\l^ ■ 

This entails {}jjhf) * (Ph'^b) f*{plB) in i° and {■ijjhf) * iPhls'^) f*{plB'^), 
and finally {tjjhf) *ph^ f *pm 1°". 

At last, we prove that ^p'j^f ^ f * p in L". We have indeed 



i'hf - ii'hf)*Ph 



< 



\i^hfi[x]h -y)- ^hfix - y)\\Ph{y)\^y 



< iPh'i-BWli sup Wi^hfi' + U) -'^Phflt'- +\\Ph'i-B4L°' sup 
0<u<h 0<u<h 

and these quantities vanish as h ^ 0. 



□ 
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